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The energy, momentum and temperature dependence of the quasiparticle local density of states 
(LDOS) of a two-dimensional ci^2_y2-wave superconductor with random disorder is investigated 
using the first-order T-matrix approximation. The results suggest that collective modes such as 
spin/charge density waves are relevant low-energy excitations of the cuprates that contribute to the 
observed LDOS modulations in recent scanning tunneling microscopy studies of Bi2Sr2CaCu2 0x- 

PACS numbers: 74.50.-|-r, 74.62.Dh, 74.72.-h 



One of the most widely debated issues in cuprate su- 
perconductivity is the possibihty of preformed Cooper 
pairsiiSiii and the origin of the pseudogap phe- 
nomenoni^ Recent experiments have demonstrated 
that the pseudogap phenomenon is unique to the 
hole-doped (p-type) cuprates and is absent above 
in electron-doped (n-type) cupratesiSiS Furthermore, 
in the quasiparticle tunneling spectra of the double- 
layer Bi2Sr2CaCu208-H5 (Bi-2212)ifl and the one-layer 
Bi2(Sr2-xLax)Cu206+5 (Bi-2201)ii systems, it is shown 
that the pseudogap can be distinguished from the super- 
conducting gap: the former evolves smoothly with in- 
creasing temperature whereas the latter vanishes at Tc. 
These phenomena suggest that the pseudogap may be 
associated with a competing order^'^'^ that coexists with 
the superconducting phase for T < Tc and persists above 
Tc until a pseudogap temperature T* . The competing 
quantum ordered phasei^*i^ can be manifested in the 
form of collective modes such as charge- and spin-density 
waves (CDW and SDW) in the superconducting state, 
as inferred from neutron scattering experiments in a va- 
riety of p-type cupratesii^*iS*i2*i2iiS However, whether 
these collective modes are closely correlated with super- 
conductivity remain controversial. Recent scanning tun- 
neling spectroscopic studies of the Fourier transformed 
(FT) quasiparticle local density of states (LDOS) of Bi- 
22x320,21^22 jjg^yg stimulated further discussions on the 
relevance of collective modesi^^i^'^i^^i^^i^'^i^^ While Bogoli- 
ubov quasiparticle interference apparently plays an im- 
portant role in the observed FT-LDOS in the supercon- 
ducting state, certain spectral details of the LDOS can- 
not be accounted for unless collective modes are consid- 
gj.gj^j23,24,25 jj-^ particular, the findings of 4 high-intensity 
Bragg peaks remaining above Tc in the FT-LDOS map of 
Bi-221322i cannot be reconciled with quasiparticles being 
the sole low-energy excitations. These new developments 
motivate us to reexamine the role of collective modes in 
cuprate superconductors by considering the energy (E), 
momentum transfer (q) and temperature (T) dependence 
of the resulting FT-LDOS modulations. 

We begin our model construction by noting that sub- 
stantial quasiparticle gap inhomogeneities are observed 
in the low-temperature tunneling spectroscopy of under- 



and optimally doped Bi-2212 single crystalS)2Si^ suggest- 
ing at least two types of spatially separated regions, one 
with sharp quasiparticle coherence peaks at smaller en- 
ergies Ad and the other with rounded hump-like features 
at larger energies A*. On the other hand, low-energy 
LDOS (for E < 0.5Ad) of Bi-2212 exhibit long-range 
spectral homogeneity. We therefore conjecture that dy- 
namic SDW or CDW coexist with cuprate superconduc- 
tivity and that they are only manifested in the quasipar- 
ticle LDOS when pinned by disorder. Thus, regions with 
rounded hump features in the quasiparticle spectra are 
manifestation of localized charge modulations due to pin- 
ning of collective modes by disorder, and the wavevector 
of the charge modulation is twice of that for the collinear 
SDW order, as proposed in Refs. 23-24 32. In con- 
trast, regions with sharp quasiparticle spectral peaks are 
representative of generic Bogoliubov quasiparticle spec- 
tra with a well-defined c?-wave pairing order parameter 
Ak ~ Ac; cos 20k, where A^ is the maximum gap value 
and is the angle between the quasiparticle wavevector 
k and the antinode direction. Our model therefore as- 
sumes 'puddles' of spatially confined 'pseudogap regions' 
with a quasiparticle scattering potential modulated at a 
periodicity of 4 lattice constants along the Cu-0 bonding 
directions, and the spatial modulations can be of either 
the 'checkerboard' patternS^^^ or 'charge nematic' with 
short-range stripesiiSiii In the limit of weak perturba- 
tions, we employ the first-order T-matrix approximation 
and consider a (400 x 400) sample area with either 24 
randomly distributed point impurities or 24 randomly 
distributed puddles of charge modulations that cover ap- 
proximately 6% of the sample area. For simplicity, we do 
not consider the effect of disorder on either suppressing 
the local pairing potential Arf(r) or altering the nearest- 
neighbor hopping coefficient (t) in the band structure of 
Bi-2212, although such effects reflect the internal struc- 
tures of charge modulations i 1 

Specifically, the Hamiltonian of the two-dimensional 
superconductor is given by 7i = Hbcs + Ti-imp, where 
'Hbcs denotes the unperturbed BCS Hamiltonian of the 
d-wave superconductor, Hbcs = J^wai'^k - M)CkcrCk(T + 
J2k ^k[Ck|cLkx +C-kiCkT], and Himp is the perturbation 
Hamiltonian associated with impurity-induced quasipar- 
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tide scattering potentiaL-SiSLSLS^ Using the T-matrix 
method, the Green's function Q associated with H is 
given hy = Oq + GoTGo, where Go is the Green's func- 
tion of Tide's and T — T-Limp / {\ — GoHimp) ■ The Hartrec is the tight-binding energy of the normal state of Bi- 



spin-independent (spin-dependent) interactions, Wk and 
are the Bogoliubov quasiparticle coefficients, -I- 
1, ul = [1 + (W£^k)] /2, = ek - ek 



perturbation potential for single scattering events in the 
diagonal part of Ti and for non-interacting identical 



point impurities at locations Vi is Va{(\} 
for non-magnetic {Vs) and magnetic [Vr 
whereas that for puddles with short stripe-like modula- 



2212 according to Norman et almA, ek = ii(cosfca; H 
cosfcj,)/2 -I- t2 cos /ca; cos fey -|- t3{cos2kx + cos2fcj,)/2 H 
'^iVs.m.e^'^''^^ t4(cos2kx cos ky + cos k^ cos2ky)/2 + is cos 2 fcj, cos 2 fey 



impuritiesSl, 



tions centering at r, ia^ 



^/3(q) 



^^gjq.r, 2sm(qy^^Rj)sm{q^,yRj) 



qy^x sm{2qx^y) 



(1) 



and that for checkerboard modulations is 
'2sm{qyRj)s'm{qxRj) 



E 



qy sm{2qx) 



[qx ^ qy) 



(2) 

Here all lengths are expressed in units of the lattice con- 
stant a, Rj is the averaged radius of the j-th puddle, 
and Vb denotes the magnitude of the scattering poten- 
tial by pinned collective modes. For simplicity, we ne- 
glect the energy dependence of Va.p^-y and assume that 
Vs , Vm and Vb are sufficiently small so that no resonance 
occurs m the FT-LDOS.26 For sufRciently large scatter- 
ing potentials, full T-matrix calculations become neces- 
sary as in Ref. ;27. However, large Vg^m would result in 
strong spectral asymmetry between positive and negative 
bias voltages^ which differs from experimental observa- 
tioni2ai2ii2a We also note that the energy dependence of 
Vjs^j reflects the spectral characteristics of the collective 
modes and their interaction with quasiparticles and im- 
purities. For instance, we expect ~ (^7^ for pinned 
SDW, where C is the impurity pinning strength and 7 
is the coupling amplitude of quasiparticles with SDW 
fluctuations 1-2*24 Empirically for nearly optimally doped 
Bi-2212, Rj ranges from 5 ~ lOi^ Here we take different 



values for Rj with a mean value {Rj 



10. 



Given the Hamiltonian and the scattering potentials 
Va.f3,.y{q), we find that for infinite quasiparticle lifetime 
and in the first-order T-matrix approximation, the FT 
of the LDOS p(r, E) that involves elastic scattering of 
quasiparticles from momentum k to k + q is: 



{lik+qWk (uk+qMk T Wk+qWk) 3 
+Uk+qt'k ("k+ql'k ± Wk+qUk) 3^ 
+Wk+qMk (Wk+qt^k ± Wk+q^k) 3 
-Wk+qWk (Wk+qMk T t^k+q^k) 3 



X 



(cj-Bk+ii5)(i^--Ek+q+»i5) 
1 



{uj+Ei,+iS}(uj-Et^+q+iS) 
1 



1 



ti_5 = -0.5951,0.1636,-0.0519,- 
is the chemical potential, and ~ 



iLj+Ek+iS){u+Ek+q+iS} 



}3) 



Here N is the total number of unit cells in the sam- 
ple, and 9[. . . ] denotes the imaginary part of the quan- 
tity within the brackets, which is related to the equal- 
energy quasiparticle joint density of states. The up- 
per (lower) sign in the coherence factor applies to 
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FIG. 1: Calculated energy-dependent Fourier transform (FT) 
maps of quasiparticle LDOS in the first Brillouin zone 
with randomly distributed non-magnetic point defects us- 
ing Eq. and 14: (a) = 40 meV and (a;/ Ad) = 
±0.15, ±0.45, ±0.75 (up and down from left to right); (b) 
Ad = 20 meV and (a;/Ad) = 0.15,0.45,0.75 (left to right), 
(c) Schematic illustration of the equal-energy contours and 
representative modulation wavevecto rs q 4 , q r and q^, which 
correspond to q^, q^ and qj in Refs. l2(12li 

Using Eq. Q and Va.fi ^-^(c^., we obtain the energy- 
dependent FT-LDOS maps in the first Brillouin zone for 
non-magnetic point impurities in Fig. ^ with two differ- 
ent Ad values and for pinned SDW (with spin-dependent 
coherence factor) in Fig. |21 whereas the corresponding 
LDOS modulations due to Va,/3^^(q) in real space is 
shown in Figs.|2Ia)-(c). For non-magnetic point impurity 
scattering at T <^Tc, the intensities associated with 
and qp are much stronger than those of q^ , as shown in 
Fig- n and also in Fig. QJa). The results in Fig. differ 
from the STM observation22iSi that reveals comparable 
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FIG. 2: Energy-dependent FT-LDOS maps with randomly 
distributed pinned SDW using Eq. ^ and V.,: (a) Ad = 40 
meV and {uj/ Ad) = ±0.15, ±0.45, ±75, up and down from left 
to right; (b) Ad = 20 meV and (u/Aa) = 0.15, 0.45, 0.75, from 
left to right. The FT-LDOS does not exhibit discernible dif- 
ferences in the spectral characteristics except the total inten- 
sities if we simply replace V-y by V/s and assume non-magnetic 
coherence factors in Eq. Q. 



intensities associated with and , and weaker inten- 
sities with q^. Interestingly, the intensities of q^ and 
Qs c become reversed if one assumes magnetic point im- 
purity scattering, as illustrated in Fig. 0Jb). However, 
there is no evidence of magnetic scattering in the samples 
used in Refs. f20 21. In contrast, the presence of pinned 
collective modes, regardless of CDW or SDW, gives rise 
to much stronger intensities for q^ (by about two orders 
of magnitude), as shown in Fig. |21 Thus, the empirical 
FT-LDOS mapsSSiSi cannot be solely attributed to quasi- 
particle scattering by non-magnetic point impurities. 
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FIG. 3: Real space quasiparticle LDOS for a (400 x 400) 
area at T = due to scattering by (a)non-magnetic point 
impurities, (b) pinned CDW and (c) pinned SDW, for Ad — 
40 meV and lu = 30 meV. 

The relevance of collective modes become indisputable 
when we consider the temperature dependence of the FT- 
LDOS. As shown in Fig. Ela), in the hmit of T -» T", 
the g- values contribute to the FT-LDOS map become sig- 
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FIG. 4: Evolution of the relative intensities of FT-LDOS with 
energy (lu) for qa, Qb and qc as defined in Fig. 1(c) and Vs, Vm 
and Vo all taken to be unity: quasiparticle scattering by (a) 
single non-magnetic point impurity, and (b) single magnetic 
point impurity. 
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FIG. 5: The FT-LDOS maps at T = 0, 0.75Tc and Tc (from 
left to right) for (a) point impurities Va{q) and (b) pinned 
SDW V.y{q). We assume Ad(r) = Ad(0)[l - {T/n)Y^\ 
Ad{0) = 40 meV, tunneling biased voltage — 18 mV, 
and Tc — 80 K. Besides temperature dependent coher- 
ence factors, the thermal smearing of quasiparticle tunnel- 
ing conductance (dl/dV) is obtained by using (dl/dV) cx 
1/ pq{E)(df /dE)\(^E-eV)dE\, where f{E) denotes the Fermi 
function, (c) I^aI-vs.-I^ (biased voltage) dispersion relation 
for pinned SDW at T = and Tc. 



nificantly extended and smeared for point-impurity scat- 
tering. In contrast, pinned SDW yields strong intensi- 
ties in the FT-LDOS map only at q^ for T >~ Tc, as 
shown in Fig. |SJb) . The overall energy dispersion due to 
pinned SDW is weaker than that due to point impuri- 
ties, as shown in Fig. [SJc) for |q^|-vs.-y (biased voltage) 
at both T = and T = Tc. In particular, we note that 
the dispersion is further reduced at Tc. These findings 
are consistent with recent experimental observation by 
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The energy, momentum and temperature dependence 
of our calculated FT-LDOS in Figs. I1I5I is supportive 
of spatially modulated collective modes being relevant 
low energy excitations in cuprates besides quasiparticles. 
In particular, only pinned collective modes can account 
for the observation in the FT-LDOS map above Tc. Al- 
though our simplified model cannot exclude CDW, we 
note that pinned CDW would have coupled directly to 
the quasiparticle spectra and resulted in stripe-like peri- 
odic local conductance modulation, which has not been 
observed in STM studies. On the other hand, various 
puzzling phenomena seem reconcilable with the SDW 
scenario. For instance, the nano-scale gap variations ob- 
served in 61-2212"^^ may be understood by noting that 
the LDOS in regions with disorder-pinned SDW contains 
information of disorder potential coupled with quasipar- 
ticles and SDW, so that the hump-likc spectral features 
at ±A* represent neither the SDW gap nor the super- 
conducting gap Ad, and the values of A* vary in accor- 
dance with the disorder potential. The long-range spatial 
homogeneity of quasiparticle spectra in YBa2Cu307_5 
(YBCO)'^^ as opposed to the strong spatial inhomogene- 
ity in Bi-2212 can also be reconciled in a similar context. 
That is, SDW can be much better pinned in extreme two- 
dimensional cuprates like Bi-2212 than in more three- 



dimensional cuprates such as YBCO. Furthermore, SDW 
can be stabilized by magnetic fields )^2*2^i^ which natu- 
rally account for the checkerboard-like spectral modula- 
tions around the vortex cores of Bi-2212. ^^''^'^ Finally, the 
smooth evolution of the pseudogap phase with tempera- 
ture through Tc may contribute to the anomalously large 
Nernst effect under a c-axis magnetic field above 
with spin fluctuations responsible for the excess entropy. 

In summary, we employ first-order T-matrix approx- 
imation to study modulations in the quasiparticle FT- 
LDOS of cuprates as a function of energy, momentum and 
temperature. Our results suggest that a full account for 
all aspects of experimental observation below Tc must in- 
clude collective modes as relevant low energy excitations 
besides quasiparticles, and that only collective modes can 
account for the observed FT-LDOS above Tc. 
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